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ABSTRACT 

I give analytic expressions for the relative number, luminosity, and mass density of 
disc galaxies as a function of surface brightness. These surface brightness distributions 
are asymmetric, with long tails to lower surface brightnesses. This asymmetry induces 
systematic errors in most determinations of the galaxy luminosity function. Galaxies 
of low surface brightness exist in large numbers, but the additional contribution to the 
integrated luminosity density is modest, probably 10 - 30%. 
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1 INTRODUCTION 



The space density of Low Surface Brightness (LSB) galaxies has long been a contro- 
versial and confusing subject. It is basic to our inventory of the contents of the universe, 
and is crucial to many aspects of extragalactic astronomy. For example, the density of 
LSB galaxies has an impact on the luminosity function, faint galaxy number counts, Lya 
absorption systems, and theories of galaxy formation. In this paper I derive analytic 
expressions which quantify the density of discs of all surface brightnesses. 

In order to clarify some of the points of confusion, it is important to quantify what 
is meant by the term 'LSB'. Different people use it to mean different things, leading to 
contradictory interpretations of otherwise consistent data. One frequent but misleading 
use is treating 'LSB' and 'dwarf as synonymous. Dwarf galaxies may be a subset of 
galaxies which are low in surface brightness, but there exist many LSB galaxies which can 
not be considered dwarfs either in terms of size or luminosity (Romanishin, Strom & Strom 
1983; Phillipps & Disney 1986; Bothun et al. 1987; Davies et al. 1988; Impey & Bothun 
1989; Bothun et al. 1990, Irwin et al. 1990; Knezek 1993; Sprayberry et al. 1993; McGaugh 
& Bothun 1994; de Blok, van der Hulst, & Bothun 1995a). Some seem to refer to LSB 
galaxies as things which are totally invisible optically. Others use the term to refer to any 
galaxy which is fainter than some canonical value, like that of Freeman (1970) or the night 
sky. 

In order to provide a consistent definition, it is necessary to quantify the luminosity 
profiles of galaxies. Sersic (1969) found that the radial luminosity distribution could be fit 
by 

E(r) = Eoe-('-/'*)' (1) 

where Eq is the central surface luminosity density (in LqPC~'^), h is a physical scale length 
(in kpc), and /3 is an index which varies the shape of the profile. This is a generalization 
of the exponential disc (/3 = 1, de Vaucoulers 1959) and r^^'^ (/3 = |, de Vaucoulers 1948) 
profiles. 

Since disc and bulge are believed to be distinct physical components, and most debate 
concerns the numbers of LSB disc galaxies, I will limit the discussion to profiles with /3 = 1 
(i.e., spiral, irregular, and dE galaxies). It is straightforward to generalise the results 
presented below for different profile shapes, but in principle a multiparamter approach 
which adequately represents the full appearance of galaxies is required (see McGaugh, 
Bothun, & Schombert 1995a). 

In magnitude units the exponential profile is 

/i(r) = /io + 1.086-, (2) 

a 

where /io is the central surface brightness and a is the angular scale length projected on 
the sky by a galaxy of physical scale length h at distance d. The surface brightness and 
size of a disc galaxy are characterised by /io and h. For pure disc systems, /io is simply 
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related to the effective surface brightness. I wiU therefore use /xq to characterise the global 
surface brightness of disc galaxies. 

Using /io it is now possible to define what is meant by low and high surface bright- 
ness (HSB). Freeman (1970) asserted that all spiral discs had the same central surface 
brightness, 

/io = 21.65 lb 0.30 B mag arcsec"^. (3) 

This has become widely known as 'Freeman's law', and is often used as a working definition 
of HSB. To be slightly more specific, I will define as high surface brightness any disc with 
jiQ < 22 B mag arcsec"^, which corresponds to Eq = 100 LqPc"^. Any disc which satisfies 
this criterion is likely also to fall within the narrow range specified by Freeman's law, and 
the particular value is chosen to bracket it on the faint side. On the bright side, I will 
describe as very high surface brightness (VHSB) anything which fails to obey Freeman's 
law by being too bright, with jiq < 21.25 mag arcsec"^ or Eq > 200 LqPc"^. 

The most obvious dividing line for defining LSB is the brightness of the night sky, one 
magnitude fainter than the Freeman value in the best of conditions. To give a linear factor 
of two between this and HSB, and to be clearly fainter than the night sky, let us describe 
as LSB any disc galaxy with jiq > 22.75 mag arcsec"^. Fig. 1 shows examples of idealised 
high and low surface brightness galaxy profiles. 

Though I would prefer to avoid the pollution of excessive nomenclature, further quan- 
titative distinctions are useful. The transition region 22 < /io < 22.75 between HSB and 
LSB I will call intermediate surface brightness (ISB), since objects in this range do not 
obey Freeman's law but are not particularly dim either. It is also useful to denote the 
extremes of the LSB realm by very low surface brightness (VLSB: 24.5 < /io < 27) and 
extremely low surface brightness (ELSB: /io > 27 mag arcsec"^). Most LSB galaxies which 
are known in the field are too bright for either of these categories, though some examples 
do exist. Nonetheless, it is possible to see VLSB galaxies on survey plates if one looks hard 
enough. It is not until the ELSB regime that galaxies become effectively invisible to the 
deepest existing photographic surveys; only one galaxy this faint has so far been reported 
(GP 1444; Davies, PhiUipps, & Disney 1988). 

The nomenclature utilised here is summarized in Table 1. While such quantitative 
definitions are necessary to make sensible use of terms like LSB, it should be realised that 
there is some distribution of central surface brightnesses. Obtaining this should be our 
goal. 

2 THE RELATIVE NUMBERS OF LSB GALAXIES 

To begin in this direction, it is possible to make a crude estimate of the relative number 
of galaxies in two widely separated surface brightness bins. Recent surveys (Schombert 
& Bothun 1988; Schombert et al. 1992) have identified large numbers of LSB galaxies by 
visual examination of the plates of the new Palomar Sky Survey. The selection criterion of 
the Schombert et al. (1992) sample (henceforth referred to as the LSB sample) mimic that 
of the diameter limited UGC (Nilson 1973): 6 > 1'. The only difference is the greater depth 
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and lower grain noise of the new plates. Objects are not included if already catalogued in 
the UGC. 

Extending the UGC in this way will incorporate galaxies which just missed the original 
diameter criterion. This will include HSB galaxies which are on average slightly more 
distant as well as lower surface brightness galaxies. In order to reject the former and 
isolate from the surface brightness distribution some portion which is truly LSB, only 
galaxies were included which were essentially undetected on the POSS-I. This guarantees 
that only shallow profile, LSB galaxies are selected. Surface photometry of a subset of 
galaxies from the LSB sample shows that this procedure selects a population which has a 
typical central surface brightness of /xq = 23.4 B mag arcsec"^ (McGaugh & Bothun 1994; 
de Blok et al. 1995a). The distribution is very sharply peaked around this value, with a 
standard deviation of 0.6 and a high kurtosis. 

Note that /xq = 23.4 is ~ 5a fainter than the Freeman value. If Freeman's law were 
correct, and the surface brightness distribution really was a sharply peaked function, then 
these LSB galaxies should not exist. Yet they represent an ~ 11% enhancement over 
the number of galaxies in the UGC (Schombert & Bothun 1988), and the UGC already 
contains a fair number of LSB galaxies (Romanishin et al. 1983; McGaugh & Bothun 1994; 
de Blok et al. 1995a). 

For a direct comparison to bona-fide HSB galaxies, a sample is required which has 
surface photometry so that the adopted definition of HSB can be quantitatively satisfied. 
Such a sample exists in the data of van der Kruit (1987), who claimed to confirm Freeman's 
law for large, early type discs. This HSB sample contains galaxies which have /xq = 
21.5 lb 0.4 i?j mag arcsec"^, the difference between this and the Freeman (1970) value 
being exactly compensated by the difference between the B and Bj bands. There are 74 
disc galaxies larger than 2' over 450 square degrees of sky, of which 51 were sufficiently face 
on to perform surface photometry. Of these, 37 qualify as Freeman discs, the remainder 
being classified as LSB dwarfs. Assuming that an equivalent fraction of the discs without 
surface photometry are HSB gives a surface density of Freeman discs of A/hsb ~ 0.12 per 
square degree. 

In comparison, Schombert et al. (1992) found 198 LSB galaxies larger than 1' over 
2375 square degrees of sky. Thirteen of these are classified as peculiar ellipticals (with 
shells or some other LSB extension), leaving 185 bona fide LSB galaxies. Only seven of 
these are dEs, the rest are spirals and irregulars. Since dEs are low surface brightness and 
generally have exponential profiles, I retain them but obviously this makes no difference. 
It is important that the sample consists mostly of spiral galaxies, as it is often stated that 
Freeman's law does not apply to dwarfs. This term is rarely quantified in this context, 
and can follow from any of morphology, size, or luminosity. The LSB sample does not 
consist of dwarfs by any of these definitions. The surface density of LSB spiral galaxies is 
A/lsb ~ 0.08 per square degree. 

An estimate of the relative space density can be made by noting that the actual 
number density n is related to the surface density by n = OA/'/U, where O is the solid 
angle and V is the volume sampled by a survey. The ratio V/fl can be estimated for 
these diameter limited surveys from equation (2) given values of the characteristic /xq of 
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each sample and the isophotal level /i^ at which the diameter 6 = 2r = 1.84a(/i^ ~ A*o) is 
measured. Since a = hid, the maximum distance at which a galaxy can lie and meet the 
survey requirement ^ > is 

h 

(/= 1.84— (/i^ -/io), (4) 
so the volume probed by the survey is 

V=^-^^d\ (5) 



This suffices for the HSB survey, but not for the LSB survey because galaxies already 
catalogued by the UGC are specifically excluded. This places a minimum as well as max- 
imum limit on the volume probed; to calculate this properly the volume already surveyed 
by the UGC must be subtracted off. In effect, the LSB survey only probes a shell of width 

^ = Jlsb -t/uGC = 1.84^(/iLSB _^yGC^ (6) 

stemming from the difference in the isophotal levels at which the diameters are measured. 
From CCD data, the diameters of the LSB survey are determined to be measured at 
= 26.0 ± 0.3 mag arcsec'^ (Schombert & Bothun 1988; McGaugh 1992), and those 
in the UGC at /xy^"^ = 25.3 ± 0.7 (CorneU et al. 1987; see also Paturel 1975; Fouque & 
Paturel 1985). The volume sampled by the LSB survey is thus 

U = 47rO / R^dR= —Q{3d'^S -idS"^ + S^) (7) 
Jd-s 3 



which reduces to the familiar approximation V = ArrQd'^S for 6<id. 

The relative number density of LSB to HSB galaxies now follows from the observed 
surface densities and the volumes sampled by each survey: nLSB/?^HSB = A/lsb/A/hsb x 
(U/0)hsb/(U/0)lsb- Substituting the above expressions for U/0 gives 

nLSB ^ ATl fe\V rhn\ 



where for brevity of superscripts the value specific to each catalog is denoted by L for 
LSB, H for HSB, and U for UGC. All of the relevant values have been discussed, except 
for the ratio of scale lengths /insB/^LSB- The LSB catalog is dominated by objects with 
scale lengths similar to the HSB galaxies studied by van der Kruit (1987), so I assume 
hnsB/hLSB = 1 (Schombert et al. 1992; McGaugh & Bothun 1994; de Blok et al. 1995a). 




5 



If instead HSB galaxies are on average larger than LSB galaxies (/iHSs/^LSB > 1), this 
will increase the relative number of LSB galaxies since nLSB/?^HSB oc (/insB/^LSB)^- 

Inserting the relevant values (A/^sB = 0.08, A^ISB = 0.12, 6^^^ = 1, ef^^ = 2, 
hRSB/hLSB = 1, l^f^^ = 21.5, //LSB = 23.4, iif^^ = 26.5, /x^sb = 26.0, /xy^^ = 25.3) into 
equation (8) yields 

log ( ^) = -0.03. (9) 

V'^HSB / 

That is, there are approximately the same number of galaxies with /xq = 23.4 as with 
lio = 21.65. 

This result is obviously very different from the notion that all spiral discs have the 
same central surface brightness. Integrating the sharply peaked Gaussian distributions of 
Freeman (1970) and van der Kruit (1987) to obtain the density expected for galaxies with 
jiQ > 23.4 mag arcsec"^ (i-e-, those > 5cr deviant). Freeman's law predicts 

log ( = -6.52. (10) 

V'^HSB / 

There are several uncertainties in this calculation, but none which come close to six and a 
half orders of magnitude. Aside from the ratio of scale lengths, variation of the parameters 
in equation (8) over reasonable ranges does not cause the answer to change by more than 
a factor of a few (see Table 2). Even if several factors conspire together, the result is still 
not much affected. 

The variations considered in Table 2 are rather extreme. It is very unlikely that 
1-1^^^' has been misestimated by a magnitude, and it does not really matter if it has 
since the volume already surveyed for LSB galaxies by the UGC is essentially zero for 
fii < 24.3. If, on the other hand, fi^^^' is even a little fainter than 25.3, the number of 
LSB galaxies rises sharply since (5 — )• as fi^^^' — > fi^^^ ■ It is not very sensible to consider 
variations in fif^^ since all that matters to the calculation is the difference between this 
and fj.^^^ . This difference is not sensitive to zero point errors, only to improperly corrected 
nonlinearities. The selection isophote of the LSB sample is well determined and uniform, 
but it is possible that the subset of galaxies for which CCD photometry is available is not 
statistically representative. Even so, varying /ig does not alter the answer much, and if 
anything it is more likely that the sample as a whole is fainter than those objects which 
have been studied in detail. If so, this again raises rather than lowers the LSB number 
density. 

The ratio of typical scale lengths can make a large difference to this calculation. As 
already noted, it goes in the wrong sense required to save Freeman's law. There is no 
evidence that the distribution of scale lengths in the LSB sample of Schombert et al. 
(1992) is distinguishable from that in the HSB sample of van der Kruit (1987). It is thus 
impossible to escape the conclusion that Freeman's law is erroneous by arguing that it 
does not apply to 'dwarfs'. 
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An estimate of the error in the relative density of LSB galaxies can be made from 
counting statistics and Table 2. The 185 galaxies in the LSB sample contribute 0.07 dex 
to the error in the relative density, while the 37 bona fide Freeman discs in the HSB sample 
contribute 0.16 dex. This results in a combined statistical error of 0.18. From Table 2, a 
reasonable estimate of the errors due to fluctuations in the various /i^ is ~ 0.2 dex, and 
an equal amount comes from possible (small) deviations of the scale length ratio from 
unity. Combining all these sources of error in quadrature gives a total error of 0.34 dex. 
Using this number and equations (9) and (10), Freeman's law is formally rejected with 
a confldence of 19a. Though this error includes some estimate of the systematic as well 
as random uncertainties, the systematic problems could of course be larger. However, it 
is impossible to consider a space density of LSB galaxies less than indicated by the raw 
surface densities. The volume correction is fairly modest (a factor of ~ 10), and the LSB 
catalog is if anything less complete than the HSB one. Hence, a very hard lower limit is 
log(nLSB /it-hsb) > —1- This is still over flve orders of magnitude larger than predicted 
by Freeman's law. 



3 AN ILLUSTRATIVE SURFACE BRIGHTNESS DISTRIBUTION 



In order to understand why studies which appear to conflrm Freeman's law are so 
misleading, let us imagine a particular form for the surface brightness distribution and ask 
what we would expect to observe. As we shall see. Freeman's law arises from the failure 
to apply volume corrections to the observed apparent distribution in order to recover the 
intrinsic distribution (see Allen & Shu 1979 as well as Disney 1976; Disney & Phillipps 
1983; McGaugh et al. 1995a). 

In order to illustrate the effect surface brightness has on the volume sampled by 
a complete survey, let us hypothesise a simple distribution of disc galaxy central surface 
brightnesses, ^(/xq). I deflne ^(/io) to be the space density of galaxies at each central surface 
brightness relative to that at some flducial value /Xq: ^(/io) = ^{f^o)/n-{fJ'o)- Motivated by 
the results of §2, let us assume that galaxies exist in equal numbers faintwards of the 
flducial value. Brightwards of this, let us assume a sharp exponential cut off. Hence 



The exponential cut off at the bright end is made to be consistent with the lack of 
galaxies observed with central surface brightness brighter than the Freeman value (Allen 
& Shu 1979). The flducial value /ig corresponds to the Freeman value. To be consistent 
with a narrow apparent distribution, the exponential cut off must be sharp, 7 ^ 0.3. 

In order to illustrate the effects of /xq on diameter and magnitude selected samples 
below, I make two further simplifying assumptions. The flrst is that /xq is not correlated 
with /i, so that on average only the surface brightness matters to the selection of galaxies at 
each /iQ. The second is that complete catalogs can indeed be constructed over sufflciently 
large volumes of space that they provide a fair representation of the universe containing 
the hypothesised ^(/xq). 




for 



for 



(11) 
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3.1 Diameter Selection 



For a complete, diameter limited survey, the volume correction associated with any 
given galaxy is obtained from V/Vmax = i^e/^Y ■ Using equation (4), 6 can be decomposed 
into /i, /io , and /i^, giving the volume sampled as a function of the galaxy parameters (/xq , h) 
and the survey parameters The volume over which a particular galaxy can exist 

and satisfy the selection criteria follows from equation (5), V oc [h(iii — fio)]^. Larger 
and higher surface brightness galaxies can be detected over larger volumes, and will be 
numerically over-represented in catalogs relative to their space density. The effect is strong 
in both /io and /i, though it is stronger in /io in the sense that no galaxy with /io > f^i will be 
included in the catalog while it is in principle possible to detect galaxies with arbitrarily 
small h. Note also that V only increases as /io becomes brighter. It never reaches a 
maximum at a preferred value of /io as suggested by Disney (1976 — see McGaugh et al. 
1995a). 

Specific cases are illustrated by Fig. 1. At a particular surface brightness, changing 
the size /i by a factor of two also changes the isophotal diameter by a factor of two. So 
the relative volume over which these galaxies could be found differs by a factor of eight, 
regardless of the isophotal level at which the diameters are measured (providing that this 
is not so faint that it is beyond the edge of the discs). At a fixed size, changing the 
surface brightness /io alters the isophotal diameter in a way that depends on the difference 
between /i^ and /io. The brighter /i^, the more serious the effect. For the illustrated case of 
shifting /io by 1.5 mag arcsec"^, the ratio of isophotal diameters are ^^4^^/^24^ — 2.77 and 
Ofi^/e^^^ = 1.52. Hence the relative volumes over which galaxies of the same size can 
be detected vary by a factor of 21.25 and 3.5 for /i^ = 24 and /i^ = 26, respectively. This 
illustrates the enormous sensitivity to /i^, and the importance of measuring it consistently 
and accurately when attempting to construct complete catalogs. Note also that it is 
essentially impossible to construct a truly complete catalog when the level /i^ is chosen (as 
is usually done) to be at the ultimate limit of the survey material, since the signal to noise 
for low surface brightness images will go to zero as /io — > /i^. 

Rather than specify a particular set of survey parameters, let us focus on the relative 
functional effects of the galaxy profile by defining a fiducial galaxy of parameters /ig, h* . 
The obvious choice for /ig is the maximum in the distribution (§3.1), and h* can be chosen 
to make this an L* galaxy (as in Fig. 1). Normalised to this fiducial galaxy, the volume 
sampled by a diameter limited survey is 



h{iii - /io) 



(12) 



Properly, this should be displayed in a three dimensional plot as a surface of V which 
declines rapidly for /io > /io and h < h* (see McGaugh et al. 1995a). Since the focus 
here is on surface brightness, I shall restrict the discussion to simple two dimensional plots 
which are projections along the surface brightness axis. 
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The apparent number of galaxies observed at each surface brightness NobsifJ-o) wiU 
just be the true distribution ^(/io) convolved with the sampling function V^(/io, h). For the 
intrinsic distribution described in §3.1, 



(fij - /ip) 



for /io > /Uq 



(13) 



and 



Nobsij^o) oc 



(fij - /ip) 



g(/*o-/to)/T 



for /io < /io 



(14) 



The result is shown in Fig. 2 for several values of /i^. 

Even though the intrinsic distribution is flat faintwards of /Aq, the apparent distribution 
is very strongly peaked at /Aq, regardless of the value of /i^. All that is gained by pushing 
fii fainter is to extend the tail of the apparent distribution (cf. Allen & Shu 1979). The 
relative volume sampled becomes very small for /io two magnitudes brighter than /i^, so 
flnding any galaxies at all with /io ~ /i^ — 2 (and such objects do exist in both the UGC 
and the list of Schombert et al. 1992) immediately indicates a large space density of such 
LSB galaxies. 

3.2 Flux Selection 

For the case of galaxy selection by apparent magnitude, the survey parameters which 
need to be specifled are the flux limit and the way in which the flux is measured. Usually, 
isophotal fluxes are used. By deflnition, these are not total fluxes, nor are they simply 
related to total fluxes except in the special case that the distribution of /io is a (5-function. 
This has an important consequence: the flux by which a galaxy is selected is not directly 
related to its luminosity as with point sources, so inverting the selection function does not 
directly yield the luminosity function. 

Nevertheless, isophotal magnitudes do have the advantage that the portion of the 
total flux which they represent is rigorously deflned (i.e., that within the isophote /i^). It 
is harder to determine what is actually being measured by other flux measurement schemes. 
In principle one needs proflle information to derive the total luminosity of a galaxy, and 
isophotal magnitudes have the virtue of being formally related to this. 

The total luminosity of a galaxy can be obtained from the proflle by integrating 
equation (1). For disc galaxies with /3 = 1, 




(15) 



For r 



max 



OO this is simply 



L 



27r/i2E 



(16) 
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The luminosity contained within some finite number of scale lengths x, 



X = r. 



max 



/h = 0.92(/i^ -/io), 



(17) 



IS 



Le = Loo[l - (1 + x)e-^]. 



(18) 



This is the true total luminosity for discs which truncate after x scale lengths. Discs are 
observed to remain exponential for at least three scale lengths, and usually more. This 
encompasses most of the light (Li > O.SLoo for x > 3), so the integration to infinity is a 
reasonable approximation of the true total luminosity. Equation (18) can also be used to 
relate the fiux measured within any given isophote to the total fiux. 

Galaxies with brighter /xq will emit a larger fraction of their total fiux above any give 
isophote than LSB galaxies, even of the same total luminosity. For the cases illustrated 
in Fig. 1, the fraction of the total fiux / = Fi/Foo measured by isophotal magnitudes are 
F24/F00 = 0.64 and F2&IF00 = 0.91 for HSB galaxies, and F24/F00 = 0.18 and F2&IF00 = 
0.74 for LSB galaxies. Hence, LSB galaxies will have their total fiuxes underestimated, 
and appear to be much fainter than HSB galaxies even if they have the same luminosity. 
For the present case, fhSB I fHSB = 0.28 at /i^ = 24, and fhSB I fHSB = 0.81 at /i^ = 26. 

This has a serious effect on the volume probed by magnitude limited surveys. Since 



3 /2 

V (X , at jii = 24, Vlsb/Vhsb = 0.15, and at /i^ = 26, Vlsb/Vhsb = 0.73. Hence 



both the number density and the luminosity of LSB galaxies will be underestimated by 
magnitude limited surveys in a way which depends sensitively on the effective isophotal 
level at which the magnitudes are measured. 

Though typical of large photographic galaxy surveys, these fxi are arbitrary. So too is 
the choice of /xq for LSB galaxies. Galaxies with fainter /xq are known to exist, including 
the extreme case of Malin 1 with /xq > 26 and L > L* . The important point is that 
the identification of any galaxy of faint /xq in a complete, magnitude limited sample im- 
mediately demands a large number density of such objects. There is nothing mysterious 
about this 'selection effect'; it stems simply from the fact that surveys sample a very much 
smaller volume of space for galaxies of low surface brightness. These galaxies will thus 
appear rare even if common, just as intrinsically faint stars are common even though the 
naked eye perceives mostly stars which are intrinsically bright. The difference is that LSB 
galaxies will appear rare, even if intrinsically luminous. 

To quantify these effects for arbitrary /xq and /i^, let us assume the same intrinsic 
density distribution as before. In general, the case of fiux selection is not as clean as 
diameter selection, for while the diameter is dominated by the disc, the fiux will also 
contain light from any bulge component which is not included in a simple exponential 
profile. However, the bulge is a small fraction of the total light for most disc systems, so 
the exponential profile remains a reasonable if imperfect approximation. 

The volume sampled goes as [1 — (l-|-a;)e ^]}^/^. Normalising, as 

before, to a fiducial galaxy of parameters (/io,/i*), and noting that Eq/Sq = lO""'"^*^'*''"'*''-' , 
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the relative volume sampled by a flux limited survey as a function of surface brightness 
and size is 



h 



) 



3 



1 - (1 + x)e 



— X 



-| 3/2 
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0.6(/to-/to) 



(19) 



1 - (1 + x*)e 



— X* 



For the assumed intrinsic distribution, the observed distribution NobsijJ'o) = ^(/Wo)^(/Wo) 
is shown in Fig. 3. 

The apparent distribution for flux selected catalogs is even sharper than for the diam- 
eter limited case. This is because of the exponential rather than power law dependence on 
/io- Unlike diameter limited catalogs, where the difference between /i^ and /xq is what mat- 
ters, it is the difference between /io and /ig that matters most to flux selection. Galaxies 
with low surface brightnesses will be very strongly selected against because of the faintness 
of /io, regardless of /i^. This only enters through x, and matters little. Thus, diameter 
limited surveys in principle provide a much better picture of the true galaxy population. 



With the results of §3, it is possible to understand why the issue of the number density 
of LSB galaxies, and the surface brightness distribution in general, has remained confused 
for so long. This results largely from unquantifled use of the term LSB, and from the 
misuse, or total lack, of volume corrections. In this section, I review earlier work and show 
that essentially all data sets provide consistent results when properly analysed. 

4.1 Optical Studies 

Discussion of the role of surface brightness in galaxy selection goes back at least to 
Zwicky (1957), Arp (1965), and de Vaucoulers (1974). The flrst quantitative statement 
about the distribution of galaxy surface brightnesses was in effect made by Freeman (1970). 
Complete samples did not exist at the time, so no correction for volume sampling was 
considered. 

Many people, especially Disney (1976), Disney & Phillipps (1983), Phillipps et al. 
(1987) and Davies (1990), pointed out that selection effects could cause Freeman's law. 
However, these arguments were based on the assumption of no correlation between lu- 
minosity and surface brightness. While galaxies cover a wide range in the fiQ-L plane, 
current data show that lower surface brightness galaxies do tend to be less luminous, just 
not smaller. In this case the qualitative effects are very different (see McGaugh et al. 
1995a for a detailed comparison). The same effect had already been noted by Allen & Shu 
(1979) who pointed out that Freeman's law is an overinterpretation of the data: all that 
can really be said is that galaxies brighter than the Freeman value are rare; the number 
fainter remained unconstrained. 

Freeman's law was apparently conflrmed by many workers (Freeman 1970; Schweizer 
1976; Thuan & Seitzer 1979; Boroson 1981; van der Kruit 1987), but only van der Kruit 
(1987) claimed to have a complete sample with surface photometry and further claimed 
that Freeman's law only applied to large, early type galaxies. He found a lack of galaxies 
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with /io ~ 22.5 Bj mag arcsec"^ and h ^ 2 kpc (for Hq = 100kms~"^Mpc~"^). The galaxies 
which were LSB were also small and late (Sd or later) type, opposite the results of Kent 
(1985). Nonetheless, van der Kruit (1987) concluded that selection effects as formulated 
by Disney & Phillipps (1983) were not occurring. Phillipps et al. (1987) state that van der 
Kruit (1987) misapplies their formalism, a point pursued by Davies et al. (1994). The latter 
claim to be able to reproduce the apparent distribution of the data of van der Kruit (1987) 
with the Disney & Phillipps (1983) visibility formalism if there is an implicit magnitude 
limit in his diameter limited catalog because of signal to noise constraints. However, if this 
reasoning were correct, van der Kruit (1987) would not have found the LSB dwarfs which 
he did. 

There is no obvious reason why van der Kruit (1987) missed large LSB galaxies. His 
study is unique in this point (Romanishin et al. 1983; Davies et al. 1988; Irwin et al. 1990; 
McGaugh & Bothun 1994; de Blok et al. 1995a; de Jong & van der Kruit 1994). If anything, 
the largest scale length galaxies tend to be low rather than high surface brightness (Bothun 
et al. 1990; Knezek 1993; Sprayberry et al. 1995a), the only consensus point being that 
the lack of large, HSB disks is not a selection effect (de Jong & van der Kruit 1994). The 
large diameter limit employed by van der Kruit (1987) can only survey a very limited 
volume of space for galaxies dimmer than the Freeman value, so perhaps this region does 
not constitute an adequate sample of the universe. 

Though the statistics are poor, the apparent distribution of central surface bright- 
nesses in the data of van der Kruit (1987) is roughly what is expected from Fig. 2. There 
is one large LSB galaxy in his sample (NGC 4392) which, when corrected for volume sam- 
pling, implies a greater density than all the large Freeman discs combined (see his Fig. 6). 
Since there was only one such galaxy in his sample, van der Kruit (1987) chose to ignore 
it. Obviously, this invalidates his primary conclusion confirming Freeman's law, even when 
restricted to large discs. 

Aside from the possibility of Freeman's law being a real effect or due to selection 
effects, it has also been suggested that it could be caused by improper disc-bulge deconvo- 
lution (Kormendy 1977; Phillipps & Disney 1983; Schombert & Bothun 1987; Davies 1990; 
Ronnback 1993) or inclination and internal extinction (e.g., Peletier & Willner 1992; see 
also Giovanelli et al. 1995). While important, neither of these effects are really pertinent 
here. They are usually small, and almost never as large as the 2 mag. difference between 
Freeman's law and the LSB galaxies of Schombert et al. (1992). Moreover, they do not 
matter much to selection, which is by far the dominant effect. 

More recently, Bosma & Freeman (1993) and Roukema & Peterson (1995) have ad- 
dressed this problem directly with further surveys. The results of Bosma & Freeman (1993) 
are examined in detail below. Roukema & Peterson (1995) claim to find a small number 
(~ 6%) of LSB galaxies. However, they fail to quantitatively define LSB, or to isolate 
a population which could be described as such (though their Fig. 5 does illustrate the 
absurdity of Freeman's law). Their selection criteria are subjective, as is their estimate 
of the redshift of completeness. This has a strong influence on the results since V 
— a no less arbitrary choice of Zc would give a number of 100%. In essence, they make 
both major mistakes which have befuddled this fleld: no quantitative deflnition of terms. 



12 



and improper correction for volume sampling. When these problems are considered, es- 
sentially all published data are consistent despite the many contradictory and misleading 
interpretations given them. 

4.2 The Diameter Distribution 

Bosma & Freeman (1993) approach the problem of the central surface brightness 
distribution by examining the distribution of diameter ratios on plates of different depth. 
They measure only the apparent distribution, without correction for the volume sampled. 
In this section, I derive this correction and apply it to their results. 

The ratio of the diameter of the same galaxy measured on two plates of different 
limiting isophotes /i^^ and fxi^ is simply related to the central surface brightness /io- In 
effect, this is a very crude form a surface photometry which fits a straight line to two 
points to estimate /io- Bosma & Freeman (1993) define the parameter F to be the observed 
diameter ratio, which is related to /io by 

F = — = ~ ^° (20) 

^£2 I^i2 ~ 1^0 



The samples are selected by diameter on the deeper plates (by choice /i^^ > f^i^)- 
Bosma & Freeman (1993) employ plates from the first Palomar sky survey ('Pal'), the 
ESO-B survey ('ESO') and SRC-J survey ('SRC'). They derive values for the limiting 
isophotes of these materials of /if"' = 24.62, fif^'^ = 25.10, and /if^*^ = 25.59 essentially 
by demanding that Freeman's law be true. This causes their derived values to differ 
significantly from the actual measured values: /if'*' = 25.3 (Cornell et al. 1987), /if"^*^ ~ 26 
(Lauberts & Valentijn 1989), and /if^*^ ~ 27 (Corwin, de Vaucoulers, & de Vaucoulers 
1980). However, the numerical value of fxi is not as important as their failure to consider 
what one should expect to see. 

The expected apparent distribution for the intrinsic model of §3.1 can be derived by 
substituting F for /io in equations (13) and (14). From the definition of F, 

1^0 = p _ ^ — • (21) 
For a flat intrinsic surface brightness distribution, the apparent distribution will be 

r.^--,.. ) . (22) 

fj'ii r*— 1 / 

where F* corresponds to /Aq. This is plotted in Fig. 4 together with the data of Bosma & 
Freeman (1993). 

Even for a flat intrinsic distribution, the expected apparent distribution iVo6s(F) is very 
strongly peaked at F*. Indeed, the expected distribution is so sharp that it is necessary 
to consider the spread introduced by random errors in F (without which values of F < 1 
are impossible). This is complicated, because F is not directly measured, but rather is 
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the ratio of two diameter measurements. This resuhs in a dependence of or on F and 
a non-Gaussian error distribution. This, combined with the smaU expected tail towards 
large F causes a fair number of points to be scattered to high F ~ 1.5. A reasonable fit 
to the data is obtained for or ~ 0.08, which is slightly smaller than the random error 
estimate (or ~ 0.12) of Bosma & Freeman (1993). If the errors were this large, a greater 
number of galaxies would be scattered to F < 1 than actually are. 

From the apparent distribution, Bosma & Freeman (1993) conclude that 80% of spirals 
obey Freeman's law. They also claim that a quarter of galaxies possess discs which truncate 
after only two scale lengths in order to reconcile the observed spike with a Freeman law. 
These are both erroneous. In fact, the observations have precisely the form which is 
expected for an intrinsicly flat distribution. Their observation of some points with F > 2 
suggest a high density of VLSB galaxies. The identiflcation of a galaxy similar to Malin 1 
in their complete sample also implies a large density of such galaxies. 

However, F provides only a very crude measure of /xq, with errors dominating the 
entire shape of the distribution. Indeed, the sharp spike does not even provide a good 
estimate of /ig, since the error in this goes as (F* — 1)""^. Since F* is always near one, small 
uncertainties in the position of F* lead to big ones in /ig. Thus all that can really be said 
from these data is that they are consistent with an approximately flat distribution. Any 
galaxies with F > 1.3 rule out a Freeman law. 

4.3 Constraints from 21 cm Surveys 

An important constraint on the population of galaxies which are difficult to detect 
optically is provided by 21 cm surveys. These depend on the H I rather than optical prop- 
erties of galaxies, so in principle suffer no bias against gas rich but low surface brightness 
galaxies. However, sensitivity limits have restricted these surveys to very small volumes. 

The 21 cm surveys fall into two categories. The flrst type are serendipitous detections 
of signal in the off beams of pointed observations of optically known galaxies (Fisher & 
TuUy 1981; Briggs 1990). These provide the apparently most impressive constraints on gas 
rich but optically faint galaxies. The second type are blind surveys, largely independent of 
known optical galaxies (Kerr & Henning 1987; Weinberg et al. 1991; Szomoru et al. 1995). 
The latter are in principle more useful, but are to date very limited in extent. 

Strictly speaking, the flrst type of survey only constrains the mass density of 'inter- 
galactic HI clouds' (Fisher & TuUy 1981), i.e., objects which are totally invisible optically. 
If a galaxy is visible at all at the serendipitous 21 cm position, regardless of its actual 
surface brightness, it does not count towards this invisible population. The eye is very 
good at detecting extended coherent structure if told where to look, so the requirement of 
invisibility is quite restrictive, probably /xq > fJ-i- 

The diameters listed in the UGC are typically measured at the jii 25.3 mag arcsec"^ 
level, but are sometimes as faint as 27 mag arcsec"^ (Cornell et al. 1987). Hence quite LSB 
galaxies are at least discernible on the POSS-I survey prints, even if their appearance 
thereon may seriously underestimate their true size. An object must be well into the 
ELSB regime to guarantee evading optical detection entirely. Thus the flrst set of surveys 
provide no real constraint at all on the density of LSB galaxies. However, data of this sort 
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would be very useful if the associated optical identification criteria were quantified, and 
the surface brightnesses of detected galaxies actually measured. 

The case of the large H I cloud in Virgo (H I 1225+0146) is a good example of the 
problems involved in using these surveys to place limits on the LSB galaxy population. 
Originally thought to be the first discovery of an HI cloud totally devoid of an optical 
counterpart (Giovanelli & Haynes 1989), HI 1225+0146 does have an associated optical 
component (McMahon et al. 1990; Impey et al. 1990). Though the HI properties of this 
galaxy are extreme, the optical component is actually quite prominent by the standards 
of Schombert et al. (1992), having a diameter 626 ~ 2', (Salzer et al. 1991). This is twice 
the limit of Schombert et al. (1992), and HI 1225+0146 would have been incorporated 
as one of the more conspicuous LSB galaxies had that survey extended far enough south. 
Clearly, claims of strict limits on the LSB galaxy population based on these sort of data are 
meaningless without rigorous quantification of both terminology and optical search limits. 

Blind surveys are in principle much more useful, but to date cover only very small vol- 
umes. This causes large uncertainties in normalization which make interpretation of their 
results difficult (see Schade & Ferguson 1994; Ferguson & McGaugh 1995). In particular, 
these surveys fail to recover the H I mass function which is known to exist in optically 
bright galaxies. Hence their ability to constrain optically dim galaxies is suspect at best. 
(Contrast the last two sentences of the abstract of Weinberg et al. 1991). 

Szomoru et al. (1995) present more complete results from the same survey as Wein- 
berg et al. (1991), including measurements of the optical surface brightnesses of 11 of the 
HI selected galaxies. Of these, three are LSB. This is a fairly high fraction, especially 
considering that the fields were not entirely randomly selected. Those centered on known 
bright galaxies are unlikely to detect LSB galaxies, which tend to be very isolated on the 
relevant scales (Bothun et al. 1993; Mo et al. 1994). Two of the three LSB galaxies were 
discovered in blank fields. 

Obviously, not a great deal can be said with such limited statistics. Nonetheless, 
the conclusion of Szomoru et al. (1995) are consistent with those presented here. Number 
density and mass density are not the same thing; LSB galaxies could be quite numerous and 
still contain only a modest amount of mass. Moreover, as noted above, the vast majority of 
LSB galaxies should be visible on sky survey prints. Thus they do not represent some kind 
of mysterious, entirely new population, but rather a familiar if neglected one consisting 
predominantly of late morphological types which span and extend a continuum of galaxy 
properties. 

5 THE DISTRIBUTION FUNCTIONS 

In this section, I derive expressions for the relative number, luminosity, and mass 
density of disc galaxies as a function of central surface brightness. 

5.1 The Number Density 

The number density of galaxies can be extracted from the apparent distribution of 
central surface brightnesses in a complete survey by inverting the procedure described 
in §3. That is, given iVo6s(/io), one can obtain the relative number density distribution 
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from ^(/io) = NobsijJ'o) ly{jJ'o)- Unfortunately, very little data exist which meet both 
requirements of completeness and surface photometry, particularly for local samples of 
field galaxies. 

One study that does is Davies (1990). These data behave precisely as expected from 
Figs. 2 and 3 (see McGaugh et al. 1995a). The distribution ^(/io) obtained from these 
data is shown in Fig. 5. 

The distribution on either side of the peak is well fit by a straight line. It is therefore 
possible to write an analytic expression that gives the relative number density of disc 
galaxies function of surface brightness: 



The peak is at /ig = 21.9Bj mag arcsec ^, and by least squares fit giving equal weight to 
diameter and flux selected samples, the slope is 



Formally, the errors on the flt are small, ±0.07 for the slope on the faint end and ±0.2 for 
the bright end. However, systematic effects plus the relatively small numbers of galaxies 
in the fainter bins makes the actual uncertainty rather larger. 

The number density derived from the data of Schombert et al. (1992) in §2 suggest a 
faint end slope closer to m = 0. There are nearly as many galaxies in this bin at /xq = 23.4 
as in the entire sample of Davies (1990), so m = is at least as signiflcant. 

5.1.1 Correlations Between Surface Brightness and Size 

The systematic effect of greatest concern is that the volume sampling function de- 
pends on size as well as surface brightness. So far, I have assumed that size and surface 
brightness are uncorrelated, so that on average variations in h cancel out and V(fio^h) 
can be approximated by V(fio). This assumption is the natural one to make given the 
functional separation of the volume sampling into size and surface brightness terms. More 
importantly, it is motivated by the data, which show no indication of a correlation between 
these two parameters (Romanishin, Strom & Strom 1983; Davies et al. 1988; Irwin et al. 
1990; McGaugh & Bothun 1994; de Blok et al. 1995a; McGaugh, Schombert, & Bothun 
1995b). However, a weak trend one way or the other is not ruled out, so it is important 
to examine the consequences of such a relationship. 

The volume sampling depends on scale length as V{h) = for both diameter and 
flux selection. Small galaxies are hard to flnd just as dim ones are. Strictly speaking, the 
slope m describes the projection of the bivariate distribution $(/io,/i) along the surface 
brightness axis. To account for possible correlations with size, we should write 



log[(^(/io)] 



m(/io - /Uq). 



(23) 




for jio > jil 
for /io < /Uq 



(24) 




(25) 
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If size and surface brightness are correlated in the sense suggested by much unwritten 
common lore, namely, that 'dwarfs' are both smaller and lower in surface brightness than 
brighter 'giant' galaxies, then both size and surface brightness discriminate against them. 
Correcting for this additional factor due to h would require even more low surface bright- 
ness galaxies, raising the distribution in Fig. 5 and causing the fit m to be less negative. 

On the other hand, if size and surface brightness are correlated in the sense that 
lower surface brightness galaxies tend to be larger, then their increased size tends to offset 
their faint surface brightness. This would push the distribution in Fig. 5 downwards, 
towards more negative values of m. However, this would imply that the lowest surface 
brightness galaxies are quite large and luminous. Local surveys with relatively bright 
limiting isophotes would give a very misleading impression of this population. 

There is some evidence pointing in both directions. The largest disc galaxies are 
inevitably those lowest in surface brightness (Impey & Bothun 1989; Knezek 1993; Spray- 
berry et al. 1995a) so that disc galaxies exist all over the surface brightness-luminosity 
diagram up to maximums in both (McGaugh 1995). However, there may be a slight ten- 
dency for the scale length distribution to become steeper with fainter surface brightness. 
This would boost both the number of low surface brightness galaxies in Fig. 5 and steepen 
the faint end slope of the luminosity function (see Bothun et al. 1991; McGaugh 1994). 
There is not a sufficiently strong effect in either direction to alter the conclusions reached 
here, though obviously these effects will modulate the precise value of m somewhat. 

5.1.2 Cosmological Dimming 

Another systematic effect of concern for the Davies (1990) sample is cosmological 
dimming. The magnitude selected sample is limited at Bj < 19.1, where the median 
redshift is 2; ~ 0.1. At this redshift, (1 -|- zY dimming amounts to ~ 0.4 mag., plus the 
A;-correction appropriate to each galaxy. Hence, it is possible that the galaxies in the 
faintest surface brightness bins are actually higher surface brightness galaxies which have 
been dimmed to appear low surface brightness. In this case, it is obviously inappropriate 
to apply a large volume correction to derive the number of galaxies in that low surface 
brightness bin. 

In order to test the potential effects on the Davies (1990) sample, I have taken the 
empirical redshift distribution for the observed magnitude range (Koo & Kron 1992) and 
attributed the maximum amount of dimming to each galaxy. Starting in the faintest bins 
and working brightwards, each galaxy is corrected by the maximum amount allowable for 
the observed redshift distribution. This is not a statistical procedure; it is the worst case 
scenario for the analysis presented here. The results of this exercise are shown in Fig. 6. 

If maximum credit is given to cosmological dimming in this fashion, something resem- 
bling a Freeman law plus a tail to lower surface brightness can be regained. However, this 
is a by choice a pathological case, the most significant aspect of which is that it is still 
impossible to recover a true Freeman law. Indeed, a universe full of Freeman discs would 
not dim to have the observed apparent distributions (McGaugh et al. 1995a). Moreover, 
surface brightness and redshift are not observed to be correlated in any published sample, 
and the distribution derived from this contrived exercise is inconsistent with the results of 
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§2. The sample of Schombert et al. (1992) gives a high density of LSB galaxies which are 
known to be local with negligible cosmological dimming. 

Indeed, one might expect the opposite effect to be at least as important. That is, from 
the nature of the volume sampling function, the bona-fide LSB galaxies which are found in 
any complete sample should on average be of lower redshift than higher surface brightness 
galaxies simply because the latter can be seen further away. Repeating the above procedure 
in reverse, this time attributing the highest redshifts to the highest surface brightness 
galaxies, we find that the distribution is not so strongly affected. This is because the bulk 
of the observed galaxies are clustered around the same surface brightness (/Uq), and most 
of the redshift distribution is filled before the low surface brightness bins are reached. This 
stretches the surface brightness distribution, causing it to become approximately flat. This 
would bring it into agreement with the density derived from the data of Schombert et al. 
(1992), so if anything, the latter case is likely to be more important. However, neither of 
these extreme cases are likely, and mostly dimming effects probably just shuffle the high 
surface brightness bins around a bit. 

The most signiflcant effect of dimming is not the redistribution of the shape of the sur- 
face brightness distribution so much as the shift caused by the median amount of dimming. 
This is ^ 0.4 mag., depending on the colours of the galaxies. This is the amount needed 
to bring the value of /Xq determined from the Davies (1990) observations into agreement 
with local samples (van der Kruit 1987). 

Hence, the surface brightness distribution ^(/io) is well described by two straight lines 
which meet at ji^ = 21.5 Bj mag arcsec"^. The slope of the faint end is probably in the 
range — 0.3 < m < with the latter value perhaps somewhat more likely. The slope is much 
steeper on the bright end, with m 2.6. (Note that this corresponds to 7 = 0.17, sharper 
even than assumed in §3.) Signiflcant uncertainties remain in all of these parameters, and 
even /Xq is not tremendously well determined. 

5.1.3 Recent Data 

Since the initial submission and presentation of these results (McGaugh 1993), much 
work has been done which strongly conflrms them. Schwartzenberg et al. (1995) flnd a 
distribution which is flat until a sudden dramatic rise in the numbers of VLSB galaxies. 
Their data are very deep images obtained with the AAT f/1 camera, and the result is 
sensitive to dimming effects. Dalcanton (1995) identifled a large number of VLSB galaxies 
in deep drift scans consistent with a flat distribution, but is obliged to make the same 
assumption about the distribution of scale lengths as I have done. Though not quite as 
deep, there are two studies (Sprayberry 1994; de Jong 1995) which do have redshifts and 
hence need make no assumptions. Both these data sets are consistent with the distribution 
derived here, though perceptible differences exist. The data of de Jong (1995) include 
extensive CCD imaging of a sample selected from the UGC and indicate a slope very 
similar to that of the Davies (1990) data. The data of Sprayberry (1994) are selected 
from plate scans by the APM, and indicate something closer to m = 0.1. Part of the 
difference may be attributed to morphological selection: de Jong (1995) concentrates on 
spiral galaxies, as classifled by Nilson (1973). It is my qualitative impression (McGaugh 
et al. 1995b) that spiral structure tends to drop in amplitude (or at least visibility) with 
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surface brightness so that Sprayberry (1994) may simply be including more smooth LSB 
disks. The data of both Sprayberry (1994) and de Jong (1995) indicate a somewhat brighter 
/Xq and steeper cut off than found here, but these parameters are correlated and depend 
on the binning and photometric system. Given the current state of the data, all results 
appear to be consistent. 

5.2 The Luminosity Density 

Once the number density of disc galaxies as a function of surface brightness is known, 
it is straightforward to obtain the luminosity density contributed by galaxies of each central 
surface brightness. Let us define the relative luminosity density J(/io) in analogy with the 
surface brightness distribution ^(/io) so that it is normalised to the absolute luminosity 
density j* produced by /ig galaxies. That is, J(/io) = Hj^'o) I jijJ't)- If the number density 
of galaxies at each /xq is known, J is simply the product of the number times the relative 
luminosity of each galaxy £(/io,/i), i.e., J(/io) = ^(/Wo)'C(/io, ^)- The relative luminosity 
follows from equation (16), giving 



£(So,/.) = ^(-) (26) 
and 



log[£(/io, /^)] = -0.4(/io - /io*) + 21og . (27) 

Combining equations (25) and (27), the relative luminosity density is 

log[J(/io)] = (m - 0.4)(/io - /io*) - log • (28) 

This is plotted in Fig. 7 for h/h* = 1. 

Basically, this just says that even if the number density of LSB galaxies is large 
(m = 0), the luminosity density produced by progressively dimmer discs declines with 
their surface brightness. Note that J(/io) does not depend strongly on the assumption 
about h since the terms that account for the number (<f> oc h~^) and luminosity {C (x. K^) 
nearly offset. The numerical effect is stronger, so if lower surface brightness galaxies are 
on average smaller, they contribute even more to the luminosity density. 

Once J(/io) is known, it is simply a matter of the definition of terms to determine 
how much luminosity density is produced by low surface brightness galaxies. If LSB is 
taken to mean anything fainter than the Freeman value, then well over half the luminosity 
comes from low surface brightness galaxies. For the opposite extreme definition of essential 
invisibility (ELSB), almost no light is produced (< 1% by extrapolation of the trend in 
Fig. 7). 

Obviously, very different interpretations will follow from the same data if the termi- 
nology is not quantitatively defined. By the definitions adopted in Table 1, most of the 
luminosity density is produced by HSB and ISB discs, while very little is produced by 
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either VLSB or VHSB discs. A small, but significant amount (10 - 30%) is produced by 
truly LSB galaxies. 

5.3 The Mass Density 

The relative mass density p{fJ,o) follows from the product of the luminosity density 
and the mass to light ratio T. Typically, this is assumed to be the same for all galax- 
ies. However, one can do better with information about the systematics of the rotation 
properties of galaxies with surface brightness (Sprayberry et al. 1995b; Zwaan et al. 1995; 
de Blok, McGaugh, & van der Hulst 1995b; McGaugh et al. 1995c; Salpeter & Hoffman 
1995). The observations require 

T oc Eo"*/^ oc ioO-2?(''o-/*S) (29) 

where g is a bandpass dependent parameter of order unity. Given this relation between 
the mass to light ratio and the surface brightness, p = T J is simply 

log[p(/io)] = log[J(/io)] + 0.2g(/io - 111). (30) 
Substituting the above results as before yields 

log[p(/io)] = [m- 0.2(2 - q)]{po - /io) " log (31) 

which is plotted in Fig. 8 for q = 1 and h/h* = 1. 

The mass density contained in LSB galaxies is an even greater proportion of the total 
than is the luminosity density. Note, however, that the observed T-Eq relation applies to 
the mass enclosed within the edges of the discs. Extrapolating this to the total mass of the 
halos associated with the galaxies would require knowledge of the total extent of the halos 
relative to the discs. The obvious assumption is that there is no systematic dependence of 
the halo extent on surface brightness (i.e., Rhalo Rdisc)- 

The relatively large mass fraction in the halos of LSB galaxies suggested by Fig. 8 
is consistent with the predictions of the structure formation model developed by Mo et 
al. (1994) to explain the spatial distribution of LSB galaxies (see their Fig. 11) and with 
the expectations of galaxy formation theory generally (Frenk et al. 1987; McGaugh 1992; 
Bothun et al. 1993; Antonuccio-Delogu 1995; Dalcanton et al. 1995). 

5.4 The Gas Density 

In addition to the relation of enclosed dynamical mass to surface brightness, there 
exists a very similar relation for neutral gas mass fraction (de Blok et al. 1995b). Lower 
surface brightness galaxies are progressively more gas rich. This may indicate an evolu- 
tionary sequence (McGaugh 1995), with galaxies increasing in surface brightness as they 
convert gas into stars. The scatter in the (Mjj j/L)-Eo relation is rather larger than that in 
T-Eo, and it is possible that LSB galaxies with low gas content have been selected against. 
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Nevertheless, the (Mhi / L)-'Eo relation represents at least the upper envelope of the space 
occupied by galaxies. In analogy with equation (31), 

log[pgasifio)] = [m- 0.2(2 - qgasMf^o - /^o) - log (J^^ . (32) 

If there remain many undetected gas poor LSB galaxies, qgas < q and qgas could conceiv- 
ably be negative. However, from present data they are only marginally distinguishable (de 
Blok et al. 1995b) with 

qgas ~ 5 ~ 1? SO Fig. 8 suffices to display PgasijJ'o) well as p{jj,q). 



6 DISCUSSION 



The relationships set out in §5 quantify an important, previously unappreciated aspect 
of the general field galaxy population. These relations are as basic to understanding 
galaxies as the luminosity and correlation functions. As such, they provide an important 
new constraint on theories of galaxy formation and evolution. They are also relevant to 
various topical problems discussed below. 

6.1 The Luminosity Function 

It has long been recognized (Zwicky 1957; de Vaucoulers 1974; Disney 1976) that 
LSB galaxies can have an impact on determinations of the galaxy luminosity function. 
Generally, these have implicitly treated the surface brightness distribution as a (5-function 
(see McGaugh 1994; Ferguson & McGaugh 1995). This is only adequate if the surface 
brightness distribution is symmetric, so it is interesting to see what effects the actual 
<f>{po) may have. 

Galaxies of the same luminosity but different surface brightness are sampled over 
different volumes (recall Fig. 1). This means that the luminosity function can not be ob- 
tained directly from the inversion of the selection function because the selection magnitude 
is not simply related to the total magnitude. Since all surveys have an effective limiting 
isophote, direct construction of $(M) from them always begs the question of how much 
light is produced by galaxies with po > pi. However, missing LSB galaxies entirely turns 
out not to be the most serious effect. Most surveys are able to at least detect LSB and 
even VLSB galaxies, and there is little luminosity density beyond that (Fig. 7). The more 
important effect comes from the systematic underestimation of the fluxes of LSB galaxies 
(this happens by definition with isophotal magnitudes — see equation 18). This effect 
can be large (~ 2 mag., McGaugh & Bothun 1994) depending on /xq and pi (McGaugh 
1994). One expects LSB galaxies to compose a very small fraction of the total sample 
(Figs. 2 and 3), so this effect can easily go unnoticed. It is quite serious though, because 
the few detected LSB galaxies carry a volume normalised weight comparable to the entire 
rest of the survey. Small errors in their luminosity are thus magnifled enormously when 
estimating $(M). This wreaks havoc on the derived Schechter (1976) parameters (Mc- 
Gaugh 1994; Ferguson & McGaugh 1995). The real situation may be even worse than 
described in these papers, because we assumed that everything had been done correctly to 
a uniform isophotal selection level. In practice this is not the case, as pi fluctuates from 
plate to plate over a range that causes signiflcant differences in the fraction of the total 
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flux actuaUy measured. The eff'ects of varying isophotal selection levels can clearly be seen 
in the combined data of Ellis et al. (1995). Surveys with brighter limiting isophotes give 
flatter slopes a and lower normalizations ^* than do deeper surveys. How much of this 
is caused by surface brightness dependent measurement effects, large scale structure, and 
evolution is difficult to say; probably all contribute. 

It is thus difficult to estimate the degree to which the integrated luminosity density is 
underestimated by current surveys. The required information is not usually published, nor 
even extracted and retained in the original work. Fairly complete information exists for 
the APM survey (Maddox et al. 1990a), for which 24.5 < /i^ < 25. For these measurement 
isophotes the flux will begin to be signiflcantly underestimated in the ISB regime. From 
Fig. 7, I estimate that this is likely to lead to an integrated luminosity density which 
is underestimated by approximately 30%. This procedure is highly uncertain, and the 
underestimate could be as much as a factor of two (Sprayberry 1994; Dalcanton 1995), or 
as little as 10%. It is certainly not an order of magnitude or zero. 

In principle, one should construct the bivariate distribution $(/io,/i) and integrate 
over this more general quantity (see Sodre & Lahav 1993; de Jong 1995). The slope of 
the faint end of the luminosity function is particularly sensitive to this procedure. For a 
flat ^(/io), a flat {a = —1.0) luminosity function (e.g., Loveday et al. 1992) requires a size 
distribution (f>{h) oc with n = —1.35 (see Bothun et al. 1991). Most estimates of the size 
distribution (van der Kruit 1987; Hudson & Lynden-Bell 1991; de Jong 1995) are closer 
to n = —2 which corresponds to a = —1.5, but considerable uncertainty exists. Moreover, 
there is little reason to expect plausible forms of the bivariate distribution to integrate to 
precisely the Schechter (1976) form, and while this might be a tolerable approximation, 
more complicated shapes seem likely (Phillipps & Driver 1995). Ellis et al. (1995) argue 
that the slope of the local luminosity function remains flat down to at least M = —16 and 
suggest that surface brightness dependent effects must act to preserve the flat shape of 
$(M). This occurs in the case of no correlation between surface brightness and luminosity 
(model A of Ferguson & McGaugh 1995), but I flnd this much less likely than the lack 
of correlation between surface brightness and size indicated directly by the data. This 
must have an effect on the faint end slope, though the upturn need not occur immediately 
faintwards of L* . 

6.2 Faint Blue Galaxies 

The large excess of blue galaxies in faint number counts (e.g., Tyson 1988) led to many 
extreme models for galaxy evolution (e.g., Broadhurst et al. 1992). I argued (McGaugh 
1994) that the correspondence of physical properties between populations of faint and low 
surface brightness galaxies suggested a common nature, removing the need for enormous 
amounts of evolution at low redshift. However, the implication that there is a one to one 
correspondence between local LSB galaxies and distant faint blue galaxies is incorrect - 
number counts do not work this way in the cosmological context. One can only discuss 
populations statistically, with the important effects being that of differential luminosity 
density measurements between surveys and potential underestimation of the slope of the 
luminosity function (Ferguson & McGaugh 1995). 
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There is now unambiguous evidence of evolution in the galaxy luminosity function 
(Ellis et al. 1995; Lilly et al. 1995a: CFRS VI). The CFRS is particularly persuasive, 
since it is deep and uniform and Lilly et al. (1995b: CFRS I) have thoroughly considered 
the relevant effects, including those discussed here. (Note that the CFRS does not imply 
a low density of LSB galaxies as might be inferred from Fig. 9 of CFRS L Since it is 
a flux selected sample, one expects LSB galaxies to constitute a small fraction of the 
total sample even though easily detected: recall Fig. 3.) With a measurement isophote of 
jii = 28 Iab mag arcsec"^, the CFRS detects nearly all the light for galaxies well into the 
VLSB regime, sufficiently far that very little residual luminosity is missed (Fig. 7). That 
the CFRS indicates a steep faint end slope for blue objects is thus no surprise. However, it 
primarily constrains the population at z > 0.3 where the observed amount of evolution is 
reasonable, if still large. The situation locally continues to pose the most serious problems. 
An upward revision of the local (f>* as suggested by several lines of evidence (see Ellis et 
al. 1995) cures many ills, but still provides no explanation for the steep galaxy counts at 
bright magnitudes (Maddox et al. 1990b) which imply a large amount of evolution very 
recently (z < 0.1), a signiflcantly bluer galaxy population than usually assumed (Koo, 
Gronwall, & Bruzual 1993; see also Ferguson & McGaugh 1995; Babul & Ferguson 1995), 
or systematic errors (Metcalfe et al. 1995). 

The effects discussed in this paper contribute to the solution of these problems, with 
an upward revision of both (f>* and a seeming likely. However, they are not of sufficient 
magnitude to be the sole solution. A factor of 2 in (f>* is required, at the upper limit of 
what can be attributed to surface brightness measurement effects. It is unclear what the 
faint end slope actually is. It seems plausible that red, early type galaxies have a flat 
luminosity function while a steeper one would be appropriate to the late types which make 
up the excess in the counts (Marzke et al. 1994; Glazebrook et al. 1995; Driver et al. 1995), 
especially as late types tend to be ISB and LSB galaxies (McGaugh et al. 1995b). If so, 
the type dependent redshift distribution will peak at lower redshifts for later types, or 
possibly even have a bimodal distribution if the classiflcation scheme also includes strange 
things at high redshift in the late type bin. Since morphology is not quantitative, one 
really needs detailed knowledge of the trivariate distribution $(/io, /i, co/our) in order to 
interpret evolutionary or cosmological effects. If the slope of this distribution becomes 
steep for blue, low surface brightness galaxies, it would cause an 'excess' which is not well 
constrained by local surveys (cf. Phillipps & Driver 1995). 

6.3 Lya Absorption Systems 

Another quantity which requires the bivariate distribution is the cross section of galax- 
ies as absorbers along the sight of QSOs. Indeed, one really needs to know the bivariate 
distribution of the gas discs. However, Fig. 8 provides a more direct way of estimating how 
much absorption is caused by discs of various surface brightnesses, as it gives an indication 
of how much gas mass is contained in each population. Though LSB galaxies are gas rich 
as individuals, the total H I density does decline with surface brightness. 

Damped Lya and Mg II absorption systems contain most of the mass seen in absorp- 
tion. From Fig. 8, it is clear that most of these would be HSB and ISB systems readily 
visible optically. A signiflcant but small fraction should be LSB galaxies, but these too 
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should be fairly easily visible optically and so one does not expect the absorbers to remain 
unidentified when closely examined. Indeed, one example of an LSB galaxy being the 
absorber has recently been noted (Steidel et al. 1994). 

There are low column density Lya absorbers without optical counterparts (Morris 
et al. 1993). These systems do not contain much mass, but do exist in profusion. By 
extrapolation of Figures 5 and 8 into the ELSE regime, one does expect there to exist 
some such invisible galaxies, though they should contain very little mass. One might expect 
these ELSE galaxies to be weakly clustered (Mo et al. 1994), consistent with observations 
of the Lya clouds (Mo & Morris 1994). However, the extrapolations involved are huge so 
the uncertainties become enormous. At most I would say that it plausible that the low 
column density absorption systems are related to a population of ELSE galaxies. It is just 
as plausible that they are bubble-like structures unrelated to individual galaxies. 

6.4 Passband Biases: A Cautionary Tale 

The data discussed here are all selected from blue sensitive photographic plates. The 
known examples of very large, massive, VLSE galaxies similar to Malin 1 tend to be rather 
more red than the normal sized LSE population (Sprayberry et al. 1995a). This raises the 
concern that substantial numbers of massive red LSE galaxies may still be missed. 

Indeed, despite the large amounts of effort that have gone into galaxy surveys, the 
situation even in the blue remains fairly abysmal. In the red, the appropriate data do not 
exist at all. In order to examine possible passband biases, consider the case of selection on 
blue sensitive plates for galaxies of the same bolometric surface brightness but a range of 
colours. In analogy with Fig. 2, imagine that equal numbers of disc galaxies exist at each 
bolometric surface brightness. The redder discs will, for the same (/Xq"'? h), appear smaller 
on blue sensitive plates than blue discs. Hence fewer red discs will be selected at any given 

This is illustrated in Fig. 9, which shows the relative number of galaxies one would 
expected to detect as a function of B — I colour for representative combinations of bolo- 
metric central surface brightness and limiting B isophote. The results depend sensitively 
on the passband and the exact shape of the spectral energy distribution, especially as /Xq"' 
approaches /i^. There is no way to know what spectral energy distributions are really 
appropriate for the putative population of red LSE galaxies. These could be faded star 
bursts (and indeed, one would expect significant numbers of such remnants from the high 
number of known short duty cycle star burst dwarfs) or something entirely different. For 
the sake of illustration, the spectral energy distributions from the models of Euzzoni (1991) 
are utilised. 

Even for HSE galaxies selected at a fairly deep /i^, the apparent number of galaxies 
declines steadily as their colours become redder. The effect becomes more severe as either 
/io°' becomes fainter or /i^ brighter. Galaxies begin to disappear entirely for only modestly 
red colours (B — / 2) for galaxies with /ig"' ~ 23, especially at the brighter limiting 
isophotes typical of large surveys. This is several magnitudes brighter than Malin 1, so it 
is worth emphasising that we would still be unaware of this enormous galaxy had it lacked 
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a prominent bulge component. All of the known examples of galaxies similar to Malin 1 
also have prominent bulges even though most LSB galaxies have little or no bulge. 

The distribution of colours for LSB galaxies in the Schombert et al. (1992) catalog 
roughly follows the appropriate line in Fig. 9. This implies a density of red LSB galaxies 
as large as that of blue LSB galaxies (i.e., a roughly flat distribution with colour). The 
data are too limited to say more because the volume probed for red LSB galaxies by B- 
band surveys is extremely small. Considering the brightness of the night sky at redder 
wavelengths, and that contrast for blue LSB galaxies is already a problem in 5, it could 
well be that only the most prominent, observationally accessible component of the galaxy 
population has been identifled. An entire universe full of dim galaxies may remain to be 
discovered. 

7 CONCLUSIONS 

I have quantifled the relative number, luminosity, and mass density of disc galaxies 
as a function of central surface brightness. These relations can be conveniently expressed 
analytically as 

log[^(/io)] = m(/io - /io) 
for the number density, 
log[J(/io)] = (m - 0.4)(/io - /io) 
for the luminosity density, and 
log[p(/io)] = (m - 0.2)(/io - /io) 

for the mass density. This last refers to both the dynamical mass enclosed within the 
optical radius, and to the H I gas mass. See equations (25), (28), (31), and (32) for further 
details. The values of the parameters /Xq and m obtained from extant data are 

/ig = 21.5 Bj mag arcsec"^, 

m = 2.6 for /xq < /"q, and 

—0.3 < m < 0.1 for /xq > /Ug. 

These relations are consistent with all relevant data, though substantial uncertainties 
remain. They represent a major shift from the paradigm of the Freeman law which suggests 
that spiral galaxies all have nearly the same central surface brightness. The value of /Xq 
corresponds to Freeman's value, but a faint end slope of m indicates roughly equal 
numbers of galaxies at each central surface brightness fainter than /ig • This conflrms the 



25 



intuition of Disney (1976) that Freeman's law is a selection effect, but one which works 
more in the fashion described by Allen & Shu (1979). 

It is very striking that the surface brightness distribution is nearly flat. It is unclear 
what signiflcance, if any, should be attached to the particular value m = 0. The sharp cutoff 
brightwards of /ig seems analogous to the turndown in the luminosity function brightwards 
of L*, but the signiflcance of the particular value of /Xq is also unclear. What is clear is 
that these numbers are fundamental to understanding spiral galaxies. 
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FIGURE CAPTIONS 



Figure 1. Idealised galaxy luminosity profiles for four exponential discs. The galaxies are char- 
acterised by their size h and central surface brightness /io- The HSB galaxies are 
Freeman discs with jiq = 21.65 B mag arcsec"^ and the LSB disks have jiq = 23.15. 
The scale lengths are chosen to illustrate the variation of isophotal quantities at fixed 
size and luminosity for different surface brightnesses. The HSB disc with h = 3 kpc 
has an integrated luminosity of 8 X 10^ Lq, roughly L* for Hq = 100 km s~"^Mpc~ . 
The LSB disc with h = 6 kpc has this same luminosity, while the LSB disc with 
h = 3 kpc has the same physical size but only one quarter the luminosity. The small 
h = 1.5 kpc HSB disc also has L = jL* = 2 x 10^ Lq. Note that when measured at a 
particular isophotal level (e.g., the horizontal lines at fxi = 24 and fxi = 26), an LSB 
galaxy will always appear smaller than an HSB galaxy of the same size, and fainter 
than one of the same luminosity. 

Figure 2. The relative number of galaxies observed as a function of central surface brightness 
for diameter limited catalogs. The solid line shows the assumed intrinsic distribution 
with equal numbers at every jiq faintwards of the value jIq 21.5 Bj mag arcsec"^. 
Brighter than /ig the intrinsic distribution is assumed to have a sharp exponential cut 
off to be consistent with observations. The constant density assumed on the faint 
side illustrates the effects of volume sampling for diameter limited catalogs selected at 
different isophotal levels /i^. The resultant apparent distributions (dashed lines) are 
shown (normalised to Nobs = 100 at /ig) for several representative fxi (labeled). These 
are typical of shallow plate material (/i^ = 24), deep plate material (/i^ = 26), and 
very deep CCD images (/i^ = 28). The apparent distribution is very sharply peaked at 
the value /ig for all values of /i^, giving the appearance that all galaxies have /io ~ /Ug 
even if equal numbers of galaxies exist at every /io > /"o • 

Figure 3. As per Fig. 2, but for catalogs limited by isophotal magnitude. The apparent dis- 
tribution of /io is even sharper than in the diameter limited case because Nobs 
;Lo-o.6(^o-/*o) rather than (/i^ - fj^o f . 

Figure 4. The distribution of diameter ratios F, which is related to the surface brightness /io. 

The histogram is the SRC/ESO data of Bosma & Freeman (1993). The heavy solid 
line shows the expected apparent distribution for a flat intrinsic distribution. This 
is very sharply peaked because of the way F is deflned. The dotted line gives the 
expected distribution after convolution with errors, indicating that these data suggest 
an intrinsically flat distribution. 

Figure 5. The surface brightness distribution derived from the data of Davies (1990). Circles 
are data selected by diameter, while triangles are magnitude limited data. Error bars 
are from counting statistics. The lines are least squares flts to the data giving equal 
weight to points from diameter and magnitude limited samples. The distribution 
declines slowly faintwards of the Freeman value, indicating a large space density of 
low surface brightness galaxies. It cuts off sharply brightwards of /ig in analogy with 
the turndown of the luminosity function at L* . 
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Figure 6. The surface brightness distribution after maximaUy extreme corrections for the effects 
of cosmological dimming. Triangles: the flux selected data shifted by the mean cor- 
rection. Dashed line: redistribution assuming all the dimmest galaxies are the highest 
redshift ones. Dot dashed line: redistribution given the opposite assumption (see 
text). These lines give the maximum envelope to which cosmological dimming can 
modulate the surface brightness distribution derived from the data of Davies (1990). 
This can not be an important effect, as it would cause inconsistency with the data of 
Schombert et al. (1992). 

Figure 7. The luminosity density produced by disc galaxies as a function of surface brightness. 

The luminosity density peaks at the Freeman value, and declines steadily towards 
fainter surface brightnesses. Substantially more luminosity is produced by discs which 
are fainter than /ig than is produced by those which are brighter. In addition to the 
data of Davies (1990), the data of Schombert et al. (1992) are plotted as the solid 
square. For this, the solid line shows the error bar due to counting statistics alone, 
while the dashed line shows the additional estimate of systematic error. 

Figure 8. The mass density contained in disc galaxies as a function of surface brightness. Points 
as per Fig. 7. A signiflcant amount of mass is contained in the low surface brightness 
tail of the distribution in excess of that in Fig. 7 because the mass to light ratio is 
observed to increase as the surface brightness decreases. Note that LSB discs with 
/io = 23 contain ^1/3 the amount of mass that is in galaxies with /xq = /Ug = 21.5. 

Figure 9. The apparent number of galaxies of the same bolometric surface brightness that will 
be observed as a function of color. Galaxies are assumed to exist with an intrinsically 
flat distribution (equal number density at each /Xq"' and color), and are selected by 
diameter on blue sensitive plates at various isophotal levels fj^f . The plot is normalised 
to Nobs = 100 for B — I = 1.3, the typical color of a i? selected LSB galaxy. Though 
i?-band surveys provide the best contrast for LSB galaxies since the sky is reasonably 
dark there, they are quite insensitive to red LSB galaxies. The exact sensitivity 
depends on the limiting isophote and the shape of the spectral energy distribution in 
a complicated way. 
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Table 1. Nomenclature 

Name /xq Sq 

B mag arcsec"^ LqPC~'^ 

VHSB < 21.25 > 200 

HSB'* 21.25 — 22.0 100 — 200 

ISB 22.0 — 22.75 50 — 100 

LSB > 22.75^ < 50 

VLSB 24.5 — 27.0 1 — 10 

ELSB^ > 27.0 < 1 

"Satisfy Freeman's Law 
^Brightness of darkest night sky 
^PracticaUy invisible 



Table 2. Variations 
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